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I. INTRODUCTION 



Recently the behavior of vector mesons in a nuclear surrounding became a lively debated 
issue. This was especially triggered by the new CERES experiments for S-Au and Pb-Au 
collisions which show an enhancement of the dilepton yield for invariant masses somewhat 
below the vacuum mass of the p meson As argued some years ago by Brown and 

Rho this enhancement might be due to the restoration of chiral symmetry 0]. In their 
approach the masses of the vector mesons scale with the quark condensate, i.e. drop with 
rising baryonic density. Thus the p peak in the dilepton spectrum ought to be shifted to 
lower invariant mass which indeed might be an explanation for the observed enhancement 
of the dilepton yield in that region . However, also other scenarios utilizing the idea of 
chiral symmetry restoration are possible which predict a rising p mass based on the effect 
that the p becomes degenerate with its chiral partner, the al meson P]. 

On the other hand, some calculations based on purely hadronic models gave rise to an 
alternative picture. Due to modifications of the pions in nuclear medium which at least 
in vacuum form the most important decay channel of the p meson the spectral function 
of the p might be drastically changed p|-pJ]]. In addition, collisions of the vector meson 



with nucleons from the Fermi sea also influence the spectral function of the former IT4-II5 



Both pion modifications and collisions with nucleons at least lead to a broadening of the 
p peak, if not to completely new structures. Hence the enhancement in the dilepton yield 
might also be explained within a purely hadronic scenario where a lot of strength is shifted 
to lower invariant mass if the p peak becomes much broader basically without changing its 
pole position |ll3| . 

Qualitatively the Brown-Rho scenario was supported by QCD sum rule analysis which 
showed a dropping of the p mass, z/the spectral function of the p is simply modeled by a 
5-function p6|-p!8|. On the other hand, the hadronic models cited above seem to indicate 



that this approximation might no longer be appropriate if the vector meson is placed in a 
medium e.g. with the saturation density of nuclear matter. A more sophisticated ansatz 
for the spectral function of the p meson including pion corrections caused by the nuclear 
surrounding was explored in []T^^. Inserting this spectral function in the hadronic part of 



the QCD sum rule basically the same dropping of the p mass as from the simple pole ansatz 
was extracted. However, in |T^, where also p-nucleon scattering was taken into account, no 



mass shift but only a large peak broadening was found. It was shown there that even then 
agreement between the left and the right hand side of the sum rule could be obtained. This 
finding clearly casts some doubt on the quantitative predictive power of the QCD sum rule 
analysis. It seems that only after one has chosen his favorite hadronic model the sum rule 
can tell whether this model is in agreement with the operator product expansion of QCD. 

The purpose of the present paper is to explore in a systematic and as far as possible model 
independent way which spectral functions for the p meson actually are compatible with 
QCD sum rules. To do so we choose a simple Breit-Wigner parametrization for the spectral 
function treating the pole position and the on-shell width as toy parameters. Varying these 
parameters we calculate the difference of the left and right hand side of the sum rule. 
According to the finding of the two groups cited above ||20|Jl^ we do not expect to find one 
single set of reasonable parameters but at least a line or a whole band of parameter pairs. Of 
course, sophisticated hadronic models yield much more complicated spectral functions than 
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the one we will use here. However, since the spectral function appears only in an integral in 
the sum rule (cf. eq. ([T^) below) we expect that the sum rule is not sensitive to the details 
of the modeling of the spectral function, but only to its gross features. Therefore we believe 
that our analysis using a simple Breit-Wigner parametrization is also applicable to more 
sophisticated hadronic models. 

In the next section we briefly recapitulate the derivation of the QCD sum rule for the p 
meson in vacuum as well as in a nuclear medium. In section Q we present our parametriza- 
tion for the spectral function. Numerical results are given in section Finally we sum- 
marize our results in [V|. 

II. OPERATOR PRODUCT EXPANSION 

We start with the time ordered current-current correlator 

n^.(g) = I ld'xe''^^{Tj^{x)UO)) (1) 
where is the electromagnetic current of the p meson, 

i/^ = \ (^7m^ - dlf,d) . (2) 

In the following we will write down all formulae for the case of finite nuclear density pn. 
The vacuum case is recovered by simply putting Pat to zero. 

At finite baryo-chemical potential Lorentz invariance is broken. From and the baryonic 
current one can construct two independent projectors L^y{q) and T^uil) which both still 
satisfy current conservation q^L^^{q) = q^T^^iuiq) = (cf. e.g. [^). Hence the correlator 
can be decomposed in the following way: 

H^, = H^T^, + H^L^, . (3) 

The scalar functions H^ and H^ in general depend on and (f where q is defined in the 
rest frame of the nuclear medium. For simplicity we restrict ourselves here to the case of 
vanishing q, i.e. where the p meson is at rest relative to the baryonic current. In this case 
H^ and H^ become equal, 

n^(g', (f = 0) = H^(g^ cf = 0) =: H(g2) . (4) 

Actually, most of the calculations concerning QCD sum rules for vector mesons in medium 
were performed for this special case, e.g. |1I6|JT7| , |20|JI^ . Only recently there appeared growing 



interest in the explicit q dependence of the p spectral function [^,^,0. The traditional 
approach to QCD sum rules for vector mesons, i.e. the pole + continuum ansatz, was recently 
generalized to finite q by Lee . A generalization of the approach presented here to finite 
q is in progress. 

For large space like four-momenta the correlator (|l|) and thus the scalar function H can be 
calculated from the operator product expansion (OPE) of QCD [^]. Utilizing a subtracted 
dispersion relation one can make contact with the time like region by (see e.g. ||14|| ) 
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n(g2) = n(o) + cg^ + ^ [ds 



Imn(s) 

s2 (s + g2 _ ^e) 



(5) 



where we have introduced g^ = — and c is a subtraction constant. Now the basic idea 
of QCD sum rules is to calculate the l.h.s. of by the OPE and the r.h.s. by a hadronic 
spectral function. 

Following P^,|2^ the OPE for the dimensionless quantity 



R{Q' 



n(g^) 



including condensates up to dimension 6 is given by 

where SC denotes the contribution from the scalar condensates, 



(6) 



(7) 



SC 




n 



112 



(8) 



^TT^ V Try v/i-/ g4 24g4 \ TT / sig^ 

The last two contributions to (|^ arise from the non scalar condensates {Sg'-jf^D^q) and 



{Sq'jfj.DjyDcrDrq) (see |1T6|J2^ for details). The medium dependence of the condensates is 
taken into account in leading order in the baryon density piv |1^ (for a review on QCD sum 
rules in nuclear matter cf. p5|). For the scalar condensates this yields 



CTiV 



(qq) = {qq)^^^ + 

^11 In 



(9) 



and 



TT 



^G' 

TT 



O (0) 



(10) 



We take the following numerical values (cf. |T^) for 

- the strong coupling constant, = 0.36 , 

- the current quark mass, m„ = 7 MeV , 



- the vacuum two quark condensate, {qq)v; 

- the vacuum gluon condensate, ( ^-G^) 

- the nucleon sigma term, ajy = 45 MeV , 

- and the nucleon mass in the chiral limit, 



= (-250MeV)3, 
1.2 • 10-2 GeV^ 



(0) 
TV 



750 MeV . 



The factors A2 and are determined from the quark parton distributions. We use the 
values 
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A2 = 0.9 and A4 = 0.12. Finally n parametrizes the deviation of the four quark 
condensate from the product of two quark condensates, i.e. from the Hartree approximation. 
This value is still not very well determined. Therefore we will choose different values for k 



to check the sensitivity of our results. Following p!6|JT^j26[] , respectively, we take 



K 



1, 2.36, 6 



(11) 
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and neglect a possible pn dependence of n as is done usually in such studies [p!6|-[T8|,pO|,[l4 . 



This is, however, a point that clearly deserves more attention since model studies |23] show 
that in the case of finite density the RPA excitations out of positive energy states counteract 
the effects of the particle-antiparticle excitations. 

It is important to note that we have not taken into account all condensates up to di- 
mension 6. Especially, we have neglected higher twist operators since they are smaller than 
the twist two contributions PB| , PU| and their determination from the experiment is not very 
accurate yet. Their influence on the OPE for the vector mesons in medium is discussed in 

The convergence of the OPE can be improved by applying a Borel transformation to (^. 
The final result is (see e.g. for details about the Borel transformation) 



1 °° 

7rM2 j^^ Imi?HAD(s) e" 



1 + 



56 



81M6 



+ 



4M4 



mNA2PN 



5 



24M6 



m^AipN 



(12) 



where M denotes the Borel mass and we have introduced the label HAD to stress that 
Imi?HAD is calculated from a hadronic model. In addition, we have absorbed the subtraction 
constant 11(0) appearing in (^ into Imi?HAD- It is given by the Landau damping contribution 
T6| presented in the following section. 



III. PARAMETRIZATION OF THE SPECTRAL FUNCTION 



Following previous works [|1^,|1^,|1^] we decompose Imi?HAD in three parts, the contribu- 
tion from the p meson, the continuum part, and the Landau damping contribution: 

Imi?HAD(s) = 7rF^e(so - s) + ^ f 1 + ^) Q{s - s,) + 5{s)^pm , ^ • (13) 



The sum rule we will examine in the following is obtained from ([121) by inserting the decom- 
position (p!3| ) and taking the continuum and Landau damping contributions to the r.h.s.. 



56 5 

rnasKiqqf - TTm^^NA^pN ■ (14) 



81M6 24M6 
If strict vector meson dominance (VMD) is applied then S would be the spectral function 



of the p meson (see e.g. pO|). Since we are examining only a parametrization of S in the 
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following we do not have to specify whether we are using VMD or not. Anyway we will refer 
to S spectral function. 

The key issue now is the choice for the spectral function S. Clearly in vacuum the spectral 
function shows a single peak structure which is, of course, at the position of the free p mass. 
In a nuclear medium we have no guideline from experiment how the spectral function might 
look like. However, examining the hadronic functions presented in the literature (e.g. [|l3l , |T4| ) 



we also find single peak structures, possibly somewhat shifted as compared to the vacuum 
case. Thus we will also model our spectral function such that it shows one peak at invariant 
mass nip with width 7. We will treat these two quantities as free parameters to figure out 
which combinations of them are actually compatible with OPE. 
We choose the following parametrization 

^^'^^n{s-m^^ + s{r{s)y' ^^^^ 

i.e. we especially neglect a possible s dependence of the mass parameter. 

Since the integral appearing in (0) is obviously sensitive to the behavior of S{s) for small 
values of s it is important to model the threshold behavior of F in a physically reasonable 
way. In vacuum the spectral function is dominated by the decay of the p into two pions. 
Taking into account the phase space as well as the derivative coupling of the p-ir-TT vertex 
we find the threshold behavior 

rUs) -^{s- 4m2)3/2 e(. - 4ml) ■ (16) 

Keeping things as simple as possible we use for the vacuum case 



1 



7 



4m! 



4m^ 

1 T- 

\ ) 



0(s - 4m^) (17) 



with the constant 7 being the on-shell width of the vector meson. For a more sophisticated 
parametrization of the spectral function in connection with QCD sum rules in vacuum we 
refer to . 



For finite nuclear density not only the decay into pions but also the scattering with 
nucleons influences the spectral function. The latter effect also contributes below the two 
pion threshold. If the Fermi motion of the nucleons is neglected the threshold for the spectral 
function of a p meson at rest is now given by the mass of one pion since the lightest pair 
of particles which can be formed in the p-nucleon collision is a nucleon and a pion. The 
threshold behavior is dominated by the lowest possible partial wave. Without any additional 
constraint from the intermediate state formed in the p-nucleon collision we assume it to be 
an s-wave state. Hence we get 

T^,^{s)^{s-mlYl'e{s-ml) (18) 

and thus 



6 



1/2 



7 



V 



71 

s 

ml 

71 

ml 



e{s-ml) 



(19) 



If one wants to study the sum rule for very low densities both thresholds (one and two pion 
masses) are relevant. In this case one has to think about a description how to smoothly 
switch on the one pion threshold with increasing density. In the following, however, we will 
concentrate on the case of saturation density of nuclear matter. Here, if the width is chosen 
large enough, there is already a lot of strength in the spectral function below the two pion 
threshold. Hence the two pion threshold does not influence the gross features of the spectral 
function for the densities we are interested in. 

We note that there are remarkable exceptions from the single peak structure by mention- 
ing two examples: First, if the spectral function of the p in medium is modeled by a pion 
loop where the pions are dressed due to medium effects a second sharp peak at low invariant 
mass shows up besides the usual broad p peak P JT0|j20[| . Second, a two pole structure with 



two maxima appears if collisions with nucleons from the Fermi sea are taken into account 
for the calculation of the p spectral function. This additional pole is basically generated by 
the coupling of the nucleon-p system to the A^(1520) resonance Such spectral func- 

tions are clearly not covered by our toy spectral function. Therefore strictly speaking, our 
parametrization is still model dependent. However we cover a large class of models with our 
ansatz since the sum rule ([1^ is not sensitive to the details of the spectral function but only 
to its gross features like the threshold behavior. 



IV. EXAMINING THE SUM RULE 

Having set up our formalism we can ask which values for mass mp and width 7 are 
compatible with OPE, i.e. which set of parameters can be inserted in (|T3p such that the sum 
rule (0) is satisfied. Of course the sum rule is not valid for arbitrary values of the Borel 
mass M. If M becomes too small the expansion in orders of breaks down. On the 

other hand for very large values of M the contribution from perturbative QCD completely 
dominates the sum rule. In this case we are no longer sensitive to the parameters we are 



interested in. In p6| , p!8| the following recipe was suggested to determine a reasonable Borel 
window: The minimal M is determined such that the terms of order o(l/M^) on the r.h.s. of 



(p!^) contribute no more than 10% to the total value of the r.h.s. For the maximal M we 
require that the continuum part is not larger than the contribution of the spectral function 
to the l.h.s. of (P), i.e. 



^Actually, the two peak structure shows up if the spectral function is calculated in lowest order 
in the nuclear density. If higher order effects are taken into account the second peak is washed out 
to a large extent (see |15| for details). 
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OO OO , , 

[ds^(l + ^) Bis - So) e~'l^'" < Ids TiF^Qiso - s) e"^/*^' . (20) 

J OTT V TT / J S 



For given values of rup, 7, F, and Sq we calculate the relative deviation of the l.h.s. from 



the r.h.s. of (|T^) and average these deviations over the Borel window defined above, schemat- 
ically 

-"^-'max 

d= J d{M^) \l-lh.s./T.h.s.\/AM^ (21) 

mm 

with 

AM^ = - M^- . f22) 

max mm V / 

Since we are mainly interested in the ranges of reasonable values for mass and width we still 
have to determine F and sq. We use the following strategy: For given values of rup, 7, and 



So we fix F by the finite energy sum rule |16 



so so 





hence 

a. 



ds Imi?HAD(s) = lds^(^l + ^] , (23) 



IsM 



(24) 



For given values of nip and 7, sq is determined by the requirement that d in (|21|) ought to 
be minimal. 

Now we require two criteria for a pair of parameters rUp and 7 to hold. First, the difference 
d as given in (|2l| ) should be reasonably small. In the following we will explore the region of 
parameter pairs which fulfill d < 0.2% and d < 1%, respectively. Second, the Borel window 
as defined in (p^) should be reasonably large. In |2^ a window of 1.4 GeV^ was found for 



the vacuum case in the traditional sum rule approach (vanishing width). As shown in 



the Borel window shrinks with increasing nuclear density. At nuclear saturation density the 
extension of the Borel window was found to be 0.9 GeV^. For vanishing width these are, 
of course, optimal values and not lower limits. In our calculations we require the window 
AM^ to be larger than 0.9 GeV^ for the vacuum case; for nuclear saturation density we take 
a value of 0.6 GeV^ as a lower limit for AM^. 

The parameter pairs which fulfill both criteria are plotted in figs. |l],H for vacuum and 
nuclear saturation density for the values for n as given in (|Tl|). To get an idea how the two 
criteria influence the final result we have disentangled them in fig. |l| for the vacuum case 
and K = 2.36. If no constraint on the width of the Borel window is imposed the left one of 
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the dashed hnes is replaced by the dotted hne. (The other hnes stay the same.) Obviously 
the window criterion gives an important constraint by cutting off low values of the mass. 
If the lower limit for AM^ is increased the corresponding line (the left dashed one) would 
move further to the right. In all the other plots of figs. |I|,^ the line which is most to the left 
is determined by the window criterion. 

For the strict requirement d < 0.2% and for k = 1 we find only a small band of parameters 
for the vacuum case and no parameters for the case of nuclear saturation density. Note, 



however, that Hatsuda et al. [|T^,|^,0 who used n = 1 took a somewhat larger value for the 



two quark condensate. This would presumably yield a picture similar to our case k = 2.36, 
at least for vacuum (fig. |I]). Actually, only for the choice k = 2.36 the experimental values for 
the vacuum case, rrip = 0.77 GeV and 7 = 0.15 GeV, can be found inside the parameter range 
which is in accordance with d < 0.2%. This is however not surprising since n = 2.36 together 
with the values for the condensates we have used here are the result of an optimization to 
fulfill the sum rule in vacuum reasonably well The authors who have used the values 
K = 1, 6 also use different values for the other condensates. Clearly, it would be of interest 
to see how a variation of the other condensate values influences our results. For simplicity 
we have restricted ourselves here to a variation of k, only. 

The values for sq, F, and the Borel window, of course, differ for each pair of parameters 
rrip and 7. To get some idea about their magnitudes we have listed in table | their values 
for some selected parameter pairs. One can clearly observe the tendency of the continuum 
threshold sq to decrease with increasing density. This is in agreement with the results from 
the traditional sum rule approach [1T6|-[T8[] and was used as an input in the analysis of [jl4 



To demonstrate that also the values for F as determined from the finite energy sum rule 
(^) are in a reasonable range we take F = 0.01 GeV^ as derived for the vacuum case with 
K = 2.36 for the physical values nip = 0.77 GeV and 7 = 0.15 GeV (cf. table |). Assuming 
vector meson dominance F is indeed given by F ^ m'^/gp ^ 0.01 GeV^ |2^ where Qp is the p- 
TT-TT coupling constant. Also F decreases with increasing density reflecting the redistribution 
of p-strength in the medium. 

The results in figs. |l|J^ are in qualitative agreement with previous work for all values of 
K. For vanishing width the mass decreases with increasing density. This is the result of the 
traditional sum rule approach. However, with increasing width the window of "good" values 
for rup is shifted to the right, i.e. to larger values. For finite density this shift becomes even 
stronger. This confirms the finding that for a small value of the width the mass drops in 
medium , while for a large value of the width the mass might stay constant . 



With rising width the window for reasonable masses considerably grows if the weaker 
criterion < 1% is applied. Therefore, the sum rule approach loses more and more of its 
predictive power. For e.g. k, = 2.36 and d < 1% we find that for small width the mass can 
be determined reasonably well within an uncertainty of about 0.1 GeV. For large width on 
the other hand the window of allowed values for the mass grows to 0.3 — 0.4 GeV. This in 
turn means that in any hadronic model that leads to a large width for the spectral function 
of the rho meson any agreement of the calculated mass with the QCD sum rule prediction is 
not a very strong statement. In any case our analysis shows that one cannot extract model 
independent informations about the properties of the p meson from the sum rules. As 
already mentioned, in all the figures the left one of the dashed lines is basically determined 
by the second criterion which requires that the Borel window is reasonably large. If this 
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criterion is removed the left one of the dashed hnes would be placed much more to the left 
(cf. fig. in for K = 2.36). Then the uncertainty in the determination of the mass for vanishing 
width would also be much larger than 0.1 GeV. 

From the figures for different values of k we find that the masses rise with growing k. At 
first sight, one would tend to exclude the 1, 6 for vacuum since the experimental 

values for mass and width are not inside the parameter range supported by the sum rule. 
This, however, does not necessarily mean that these values for k ought to be discarded for 
finite density, too. Since k might explicitly depend on the density the values we have used 
are not unrealistic a priori. Fig. |^ shows that even for vanishing width a k which grows with 
the density can compensate for the negative mass shift otherwise caused by the density and 
restore the vacuum value of the mass. We conclude that the sum rule analysis is still plagued 
by a large uncertainty in the determination of the condensate values, especially concerning 
the four quark condensate. 



V. SUMMARY 

We have presented a systematic study of the QCD sum rule for the p meson to shed 
some light on the question which hadronic spectral functions for the p meson are compatible 
with the operator product expansion of QCD. Using a simple Breit-Wigner ansatz with 
the correct threshold behavior incorporated we explored the range of parameters for the 
p mass and on-shell width which satisfy the QCD sum rule. From the qualitative point 
of view the important message is that there is not only one pair of reasonable parameters. 
Quantitatively we find that for small (large) values of the mass also the width has to be small 
(large). Thus the lowest mass is given by the traditional calculations using a 5-function for 
the spectral function [|16|,0,p6|,IT8[ . However, especially for large values of the width a large 



range of values for the mass would also be allowed. We conclude that specific properties 
of the p meson in a nuclear medium cannot be determined from QCD sum rules without 
specifying a hadronic model. Instead the sum rules only provide a (wide) constraint for 
reasonable hadronic models. Even then, the factorization of the four-quark condensate, 
i.e. the dependence of any result on the parameter k and its possible density dependence, 
introduces an additional uncertainty. 

Nevertheless the sum rules are useful, first, to provide a constraint for reasonable hadronic 
models or, second, to incorporate additional medium dependences of hadronic parameters 
(see e.g. [^). We stress again, however, that the sum rule loses its predictive power with 
increasing width since the window for values of the mass supported by the sum rule dras- 
tically grows. For any further study concerning the comparison of a hadronic model with 
the OPE we suggest to quantify any statement about the agreement with the QCD sum 
rule. We advocate to use the criterion for a reasonably large Borel window as presented by 



Leinweber et al. p6|JT8[] together with the deviation d presented here as a tool to check the 
compliance of any hadronic model with the OPE. 

The work presented here is restricted to p mesons which are at rest in the nuclear medium. 
A generalization to finite momentum is in progress. 
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TABLES 





rup [GeV] 


7 [GeV] 


so [GeV2] 


F [GeV^] 


M2.^ [GeV2] 


[GeV2] 


d 


vacuum, k = 2.36 


opt. 


0.71 


0.01 


1.1 


0.008 


0.69 


1.6 


0.0011 


exp. 


0.77 


0.15 


1.2 


0.01 


0.69 


1.8 


0.0015 


vacuum, k = 1 


opt. 


0.71 


0.48 


0.96 


0.01 


0.51 


1.4 


0.0013 


opt.O 


0.64 


0.01 


0.95 


0.06 


0.51 


1.4 


0.0070 


exp. 


0.77 


0.15 


1.3 


0.01 


0.51 


1.9 


0.014 


vacmuii, K = 6 


opt. 


0.90 


0.08 


1.7 


0.02 


0.95 


2.5 


0.0013 


opt.O 


0.88 


0.01 


1.6 


0.02 


0.95 


2.4 


0.0013 


exp. 


0.77 


0.15 


1.1 


0.01 


0.95 


1.5 


0.0097 


nuclear saturation density, k = 2.36 


opt. 


0.66 


0.22 


0.88 


0.006 


0.54 


1.3 


0.00085 


opt.O 


0.57 


0.01 


0.78 


0.003 


0.54 


1.1 


0.0019 


nuclear saturation density, k = 1 


opt. 


0.62 


0.48 


0.68 


0.005 


0.41 


1.0 


0.0035 


opt.O 


0.49 


0.01 


0.67 


0.002 


0.41 


1.0 


0.0092 


nuclear saturation density, n = Q 


opt. 


0.74 


0.01 


1.2 


0.009 


0.74 


1.7 


0.0011 



TABLE I. Values for sq, F, M^jj^, M^^^, and d for some selected parameter pairs rUp and 7. 
opt. denotes the optimal value (minimal d) of the whole parameter range, opt. denotes the optimal 
value for the smallest width we have studied, exp. denotes the experimental values nip = 0.77 GeV 
and 7 = 0.15 GeV. 
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FIG. 1. The width 7 over the mass nip for vacuum and for different values of k. The full lines 
border the region of QCD sum rule allowed parameter pairs with d < 0.2% and AM^ > 0.9 GeV^, 
the dashed lines border the allowed region for d < 1% (same AM^). The diamond marks mass and 
width of the free p meson. The leftmost (dotted) line in the middle picture (k = 2.36) gives the 
lower boundary for the d < 1% criterion if no constraint on the Borel mass window is imposed. 
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FIG. 2. The width 7 over the mass nip for nuclear saturation density po and for different values 
of K. The full lines border the region of QCD sum rule allowed parameter pairs with d < 0.2% 
and AM^ > 0.6 GeV^, the dashed lines border the allowed region for d < 1% (same AM^). The 
diamond marks mass and width of the free p meson. 
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